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The Bowler-Sinyukov method to eliminate Coulomb interaction from a 
two-particle correlation function is discussed and tested. 
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1. Introduction 

The correlation functions of two particles with 'small' relative momenta 
provide information about space-time characteristics of particle's sources in 
high-energy nucleus-nucleus collisions [HE1E]. Within the standard 'fem- 
toscopy' method, one obtains parameters of a particle's source fitting ex- 
perimental correlation functions with theoretical ones calculated in a given 
model. Since we usually deal with electrically charged particles, observed 
two-particle correlations are strongly influenced by the Coulomb interac- 
tion. The effect of the Coulomb force is eliminated from experimental data 
by means of the so-called Bowler-Sinyukov procedure [HE]. 

The femtoscopy was applied to a large volume of experimental data on 
nucleus-nucleus collisions at SPS energy [2]. The spatial size of particle's 
sources appeared to be comparable to the expected size a fireball created 
in nucleus-nucleus collisions while the emission time of particles was signif- 
icantly shorter. It was predicted that at RHIC energies the emission time 
would be significantly longer due to the long lasting hydrodynamic evolution 
of the system created at the early stage of nucleus- nucleus collisions [HE]. 
To a big surprise the experimental data obtained at RHIC [SJ E] show a 
very little, if any, change of the space-time characteristics of a fireball when 
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compared to the SPS data. And in contradiction to hydrodynamic models 
the emission time of particles appeared to be as short as 1 fm/c. Because of 
this surprising result, which is known as the 'HBT Puzzle' |10t 111], a relia- 
bility of the femtoscopy method was questioned. Very recently it has been 
shown that the hydrodynamic calculations can be modified to give rather 
short emission time of produced particles \12\ I13j. and thus the 'HBT Puz- 
zle' seems to be resolved. Nevertheless it is still of interest to quantitatively 
check the femtoscopy method. 

Our aim here is to test the Bowler-Sinyukov correction procedure which 
is used to eliminate the Coulomb interaction from the experimental data. 
The procedure assumes that the Coulomb effects can be factorized out. The 
correction's factor is calculated for a particle's source which is spherically 
symmetric and has zero lifetime. We examine the procedure applying it to 
the computed Coulomb correlation functions of identical pions coming from 
anisotropic sources of finite lifetime. The effect of halo [T7] is also studied. 
We treat the computed Coulomb correlation functions as experimentalists 
deal with the measured correlation functions. Thus, we extract the cor- 
relation function which is supposed to be free of the Coulomb interaction. 
However, in contrast to the situation of experimentalists we know actual pa- 
rameters of particle sources which can be compared to the extracted ones. 
Our study is somewhat similar to that presented in [14J. 

We use the natural units, where c = fi = 1, and our metric convention 



We compute the correlation function using the well known Koonin for- 
mula [16J. Since the two particles of interest are described by means of 
nonrelativistic wave function, the computation is performed in the center- 
of-mass frame of the pair, as the pair motion can be treated as nonrelativistic 
in this frame. However, the source function, which gives a probability to 
emit two particles at a given space-time distance (t,r), has to be trans- 
formed to the pair center-of-mass frame (where quantities are labeled with 
asterisks). The correlation function thus equals 



where ^q»(r*) is the non-relativistic wave function of relative motion and 
£) r (r#) is the effective source function 
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where D r (i*,r*) is the 'relative' source function. Obviously, the velocity of 
the pair in its center-of-mass frame vanishes (v* = 0). The 'relative' source 
function is defined through the single-particle source function as 

D r (r,t) = Jd 3 RdTD(R-±r,T-h)D(R+±r,T + h) . (3) 

As a probability density, the source function is normalized to unity 

J d 3 rdtD(t,r) = J d 3 rdtD r (t,r) = J d 3 r D r (r) = 1 . (4) 

The Coulomb function of two non-identical particles interacting due to 
repulsive Coulomb force is well-known [15] to be 

p q (r) = e - S? T(l + i 1 !) e^ r F{ - *2 1, i(qr - qr)) , (5) 

where q = |q| and 1/rj is the Bohr radius of two-particle system which 
equals J?" 1 = 388 fm for tt-k; F denotes the hypergeometric confluent func- 
tion. As we deal with pairs of identical bosons, the wave function </? q (r) is 
symmetrized. 

We choose the gaussian form of the single-particle source function D(t, r) 
but in order to easily transform it from the source rest frame to the center- 
of-mass frame of the pair, we write it down in the Lorentz covariant form 

^ . , V detA , 1 , ,,,, . . , 

D(x) = -^r- exp[--^A^^], (6) 

where = (t, r) is the position four-vector and A^" is the Lorentz tensor 
depending on the parameters r, R x , R y and R z , which characterize the life- 
time and sizes of the source. In the source rest frame the matrix is diagonal 
with the r~ 2 , R~ 2 , R~ 2 and R~ 2 along the diagonal. The source function 
as written in Eq. @ obeys the normalization condition Q not only for the 
diagonal matrix A but for non-diagonal as well. The source function ^ is 
evidently the Lorentz scalar that is 



^V) = exp [~x'^x' v ] = exp [-\x^x u ] = D(x) 

where x'^ = L"x v and A ,fJiL/ = L fl .A ap L p u with being the matrix of 
Lorentz transformation. We note that detA' = detL detA detL -1 = detA. 

The correlation function of two identical noninteracting bosons should 
equal 2 at vanishing momentum (Cf ree (q = 0) = 2) but free correlation 
functions extracted from experimentally measured ones appear to be sig- 
nificantly smaller than 2 at q = 0. There was introduced the idea of halo 
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|17j to explain this fact. It is assumed that only a fraction / of particles 
contributing to the correlation function comes from the fireball while the re- 
maining fraction (1 — /) originates from long living resonances (0 < / < 1). 
Then, we have two sources of the particles: the small one - the fireball and 
the big one corresponding to the long living resonances. The single-particle 
source function thus equals 



where Df(t,r) and Dh(t,r) represent the fireball and halo, respectively. If 
the halo radius is so large that R^ 1 is below an experimental resolution of 
the relative momentum q, the particles coming from halo do not contribute 
to the measured correlation function and one claims that Cf ree (q = 0) = 
1 + A, where A = f 2 < 1. 

We have computed the Coulomb correlation functions for anisotropic 
gaussian sources of finite emission time. The halo has been also included. 
We use the Bertsch-Pratt coordinates [TS1 [TO] out, side, long. These are the 
Cartesian coordinates, where the direction long is chosen along the beam 
axis (z), the out is parallel to the component of the pair momentum which 
is transverse to the beam. The last direction - side - is along the vector 
product of the out and long versors. So, the vector q is decomposed into 
the q a , q s , and qi components. If the particle's velocity is chosen along the 
axis x, the out direction coincides with the direction x, the side direction 
with y and the long direction with z. 



The Coulomb effect is usually subtracted from the experimentally mea- 
sured correlation functions by means of the Bowler-Sinyukov procedure. In 
the absence of halo, procedure assumes that the correlation function can be 
expressed as 



where Cf ree (q) is the free correlation function and K{q) is the correction 
factor which depends only on q = |q|. The correction factor can be treated 
as the Coulomb correlation function of two nonidentical particles of equal 
masses and charges. The function is, however, rather unphysical as the 
pair velocity with respect to the source is assumed to vanish even so the 
calculation is performed in the rest frame of the source where the source is 
assumed to be symmetric and of zero lifetime. The correction factor K(q), 
which is described in detail in the Appendix to the paper [H] . is computed 
as 
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Fig. 1. The free correlation function Cf ICC (q , 0, 0) extracted from the Coulomb 
correlation function by means of the dilution (left) and Bowler-Sinyukov (right) 
procedure for various A. The expected free correlation function is also shown. 



where G(q) is the so-called Gamov factor equal 

G(q) = ^ ^ (10) 

W q exp(^)-l 1 * 

and D r (r) describes the spherically symmetric gaussian source of zero life- 
time and of the 'effective' radius R = \J (R^ + R'j + Ri)/3 where R a , R s and 
Ri are the actual source radii. 

To check the validity of Eq. (j5J) , we have divided the computed Coulomb 
correlation function by the Correction factor K(q). For the case of pion-pion 
correlations, the extracted free correlation function is almost identical with 
the actual correlation function of non-interaction particles. The procedure 
works very well even for strongly anisotropic sources. 

The situation is more complex when the halo is taken into account. We 
test two versions of the Bowler-Sinyukov procedure: the dilution method 
and the proper Bowler-Sinyukov one. The experimentally measured corre- 
lation functions C(q) are fitted as 



C(q) 



l-X + XK(q)) 1 + A(C free (q) - 1) 



for dilution , 

1 — A + \K(q)Cf ICC (q) for Bowler — Sinyukov . 



The correlation function Cf ree (q) extracted by means of the dilution and 
Bowler-Sinyukov procedures are shown in Fig. 1. The expected free function 
is shown for comparison. The source parameters are given in the figures. 
The parameter A is assumed to be known when C(q) are fitted. We show 
here only the function Cf ree (q) for q = (g G , 0,0) which is crucial for the 
emission time determination. As seen, the extracted correlation function is 
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distorted at small relative momenta but the width of the correlation function 
is unaltered and so are the source parameters. In our paper [20] we present 
a very detailed analysis of the Bowler-Sinyukov procedure. In particular, we 
show there that for the kaon-kaon correlations it works significantly worse 
than for the pion-pion ones. 

We conclude our study as follows. In the absence of halo the Bowler- 
Sinyukov procedure works very well for irir correlations. When the halo is 
taken into account the extracted correlation functions are distorted at small 
relative momenta but the source parameters are still reproduced accurately. 

We are very grateful to W. Broniowski, W. Florkowski, A. Kisiel and 
R. Lednicky for numerous discussions. This work was partially supported 
by Polish Ministry of Science and Higher Education under grant N202 080 
32/1843. 
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